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Abstract 

n 

Using e expansion technique proposed in we derive an effective Lagrangian (Ginzburg-Landau- 
like functional) of the degenerate unitary Fermi gas to the next-to-leading (NLO) order in e. It is 
demonstrated that for many realistic situations it is sufficient to retain leading order (LO) terms 
in the derivative expansion. The functional is used to study vortex structure in the symmetric gas, 
and interface between normal and superfluid phases in the polarized gas. The resulting surface 
free energy is about four times larger than the value previously quoted in the literature. 



1 



I. INTRODUCTION 



There has been considerable effort, both theoretical and experimental, to understand 
properties of dilute strongly interacting cold fermion systems realized, for example, in the 
experiments on trapped cold atomic gases, and in dilute neutron matter encountered in the 
neutron star crust 

mmm 

. The case of infinite scattering length and zero effective 
range of the two body interactions has been named " unitary Fermi gas" and is particularly 
interesting due to its universal properties. 



Inspired by work of Nussinov and Nussinov 



9|, an analytical technique similar to the e 



expansion in the theory of critical phenomena has recently been proposed Comparison 
with the results of the state-of-the-art Monte Carlo simulations and with experimental results 

eading order (NLO) the e expansion might be a useful 



3, 



ll| . While subsequent investigation raised doubts 



suggested that already at the next-to- 
tool in the description of this system 
about the convergence of the series [12] , one obvious way to check usefulness of the technique 
is to make predictions for various observable quantities at NLO. If the trend suggested by 
the initial NLO results of Q, 11] were to hold, it would serve as a strong encouragement 
for further investigation. 

So, as a step in this direction, in this article we derive an effective Lagrangian of the 
degenerate unitary Fermi gas in derivative expansion, and to the next-to-leading order (NLO) 
in e which is the subject of Section 2. The functional is then applied to two spatially 
inhomogeneous phenomena in the unitary Fermi gas: single vortex structure and interface 
between normal and superfiuid phases of the imbalanced gas at the critical polarization 
(Sections 3 and 4). In this paper predictions for d = 3 will be made by just setting e = 1 in 
the end of NLO calculations. More sophisticated extrapolations to d = 3 [l^ are left to 
future work. Some further applications the effective Lagrangian functional are discussed in 
Section 5. 



II. EFFECTIVE LAGRANGIAN TO NLO IN e 



The dynamics of the system is governed by the Hamiltonian 



H = ^ fylf-^- fic^j ^jja- Co^|J\^p^^|Jl^|Jl, 



(1) 
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the corresponding Lagrangian density is given by 

\ 2m J 

+ Co V'tVrV'lV'i (2) 

where ip = {'^]i'^iY'i is spin-1/2 fermion field, n = diag(/i|, yUj) are the chemical potentials 
for the two pairing species. 

The coupling cq is to be tuned to reproduce a desirable two body scattering length, a. 
In this article we will concentrate on the unitary regime, a = oo, which in the dimensional 
regularization corresponds to I/cq = [1] . It is possible to use e expansion to describe system 
near unitarity 13| , and the results of this paper may be extended to the near- unitary regime. 
This is left for future work. 

A convenient tool that we will use to study inhomogeneous configurations of the unitary 



Fermi gas is the effective Lagrangian of the theory ([2]) (see, for example, Ch. 11 of 14 1 
for an introduction). The zero temperature effective action functional r[$(x), /i, of the 
order parameter oc {ip{x)iil!{x)i), where (...) means expectation value in the ground 
state, is defined by the following 

T[^{x), n,Sij] = -^l[J{x), n,5ij] - (^j J*{x)^{x) + h.c^ , x = {t,x), (3) 

with the source J(x) satisfying 
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6J*(x 



-Q[J{x), fi,Sfi\ = -$(x) 



and 



e '^W= /D0D0Met| ^^* + S + /^T | ^ ( / (4) 



X " 



is the auxiliary di-fermion field introduced by the Hubbard-Stratonovich transfor- 



where 

mation [l|,ll3|]; I/cq = has been used; we have introduced /i = (/i|+/ij)/2, 6fi = {jj,]—jj,i)/2. 
The theory is defined in d = 4 — e spatial dimensions and the calculations are performed 
using perturbation theory in e about (i = 4 [l|, Q • 

The objective is to integrate out in (jlj) to a given order in e in the presence of J{x) 
and then perform the Legendre transformation. To approximate r[$(x), /i, which is, in 
general, non-local, we will attempt an expansion in derivatives of $(x), and, so, will assume 

r[$(x), /i, = / Lejj[$(x), /X, = / —Ve//($(x), /i, 5/i) + derivative terms. (5) 

J X J X 



So, in the case of a homogeneous medium effective Lagrangian reduces to effective potential, 
Le//[$, /i, = — Vefr($,/i) for $ = const. The effective potential for symmetric (near) 
unitary Fermi gas {6jj, = 0) has been calculated to NLO in e in [l3|; the calculation has 
been extended to NNLO in e by Arnold, Drut and Son [12^. The NLO effective potential 



for the polarized case {6^ ^ 0) was derived in [13|, ll5|. As already stated, in this article 
we will be working to NLO in e at unitarity. Within derivative expansion approach the 
task is to augment the potential with derivative terms. In the following we are going to 
argue that at least the leading order (LO) derivative terms may be determined by imposing 
gauge invariance on the theory, which allows one to avoid performing explicit Legendre 
transformation. The derivative expansion approach is reminiscent of the Ginzburg-Landau 
functional construction, but here we do not assume being in the vicinity of a second order 
phase transition. Also here we have an (approximate) control over the microscopic theory 
and are able to calculate the coefficients in the functional. Note that if one were to fail to 
find a sensible derivative expansion where only few terms have to be kept to achieve given 

ocal and a more general technique 



161]) would have to be employed. 



accuracy, then the effective action would indeed be non- 
(such as the standard Bogoliubov-de Gennes method 
But in the following we are going to argue that for a few realistic situations the derivative 
expansion is related to the e expansion and, so, a finite number of derivative terms has to 
be kept to a given order in e. 

Once the functional Q is known, the ground state is characterized by the equation of 
motion 

^^Le//[$(a;),/i,5/x] = 0. (6) 
Then, for example, free energy and particle numbers are given by 



F{fx,5fx) = -J d'^xLe//[<l(a;),/i,5/i] 



N, + N, = -^%M N,-N, = -^%M, (7) 
where ^{x) is a solution of ([6]) and Ni, N2 are the particle numbers of the two species. 
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A. Effective Lagrangian for symmetric unitary Fermi gas {5fi = 0) to NLO in e 



Let us begin by considering the case of unpolarized unitary Fermi gas {S^ = 0). The 
effective potential (that is, non- derivative terms of the effective action) for unitary Fermi 
gas has been calculated to 0{e) by Nishida and Son [l| 



Ve//($(x),/i) 



m |$(x) 

m |$(x) 
2^ 



d/2 



d/2 



\^(x) 



6 



(8) 



where 7 ~ 0.57722, is the Euler-Mascheroni constant and C ~ 0.14424. Minimization yields 



2/^o(l + £(3C-l + Log2)), /io = ^~l. 



(9) 



the NLO homogeneous medium superfluid order parameter. 



1. The Leading Order Derivative Terms 

General coordinate and conformal invariance has been used by Son and Wingate to 
constrain low energy effective Lagrangian of the unitary Fermi gas for the phonon field (the 



phase of di-fermion condensate) [iTj. However, in the situations that will be dealt with in 
this paper it will be more convenient to have an effective Lagrangian depending on both 
phase and magnitude of the order parameter. 

We will argue that within the framework of e expansion just gauge invariance will be 
sufficient to determine at least the leading order (LO) derivative terms. The basic idea is 
to consider a generalization of the theory ([2]) by gauging the U{1) particle symmetry and 
imposing gauge invariance. If a certain relationship is true for the effective Lagrangian of 
the generalized theory, then it also has to hold for the original theory. 

Our theory ([2]) is invariant under U{1)^ x U{l)i = U{l)v x U{1)a which are particle 
number symmetries for the spin up and down species. To obtain LO derivative terms for 
both polarized and unpolarized cases it is sufficient to gauge U{l)v and require that the 
Lagrangian be invariant in the presence of an arbitrary external gauge field = {Aq,A). 
Then Lagrangian ([2]) becomes 

UAoix), A(x)) = (i dt-Ao)ij- — {Vtp^ - I A4j^) -{V^ + i AtP) 

2 m 

+ CoV'JV'T^I^i- (10) 



Under an infinitesimal U{1) transformation, exp ia{x), the fields transform as 

■0^(1 + ia{x))tlj, 0^(1 + 2ia{x))(f) 

Af)-^ Ao-dta{x), A-^ A-Va{x). (11) 

For example, the case of a gas trapped in a potential U(x) is represented by Aq = — yU + U(x) 
and A = 0. To determine the LO derivative terms it will be sufficient to consider effective 
potential in the presence of = const. 

To determine the leading order (LO) derivative terms we observe that 

1. the A"^ term appears in the bare Lagrangian in the sum with Ao as 

-^UAo + ^]^; (12) 

2. the A ■ i ip'^'W ip / 2 m + h.c. term does not contribute to the coefficient of the A^ term 
in the effective action, i.e. to the Meissner mass of the U{1) gauge boson; 

3. by gauge and rotational invariance the lowest order effective action derivative term 
has the form 

Zi{<^)<^*{idt-2AQ)<^ - Z2($)(V<I>* -2^1$*) • (V$ + 2zi<l>) = 

= Zi($) <^*{idt - 2 Ao) $ - Z2($)(|V<I>P + AA^ |$n + (15) 

where ^i($) and Z2{^) are some functions. 

Functions ^i($) and Z2{^) are to be determined from the coefficient at Aq + in 
the effective potential in the background of constant gauge field A^ = {Aq,A). But the 
calculation for Aq = —fi has already been performed by Nishida and Son with the result 
shown in (JHl). Then the LO derivative terms are 

o ^ 1 ,^^,2\ 1 f l-2(7-ln2) \ fmrnV^^ 
V 4m' ' y 2|$|2e V 4 ) \ 2-k ) ^ ' 

Note that the coefficient at [A^Y term in Ve//(Ao,A) is given by a convergent integral 
and is 0{1). This will be used in the next section. 



6 



2. The Higher Order Derivative Terms 



Now let us consider the higher derivative terms. Gauge and rotational symmetries are 
not sufficient to specify derivative terms unambiguously, but as will become clear from 
the following, it will suffice to determine the order of the derivative coefficients in e. For 
simplicity, we will concentrate on time independent configurations, but a similar argument 
may be made for the terms with time derivatives. 

An n-derivative term has the form 

with various combinations of r and s > 2 compatible with the rotational invariance and such 
that 2r ■ s = n, where Zgr {j^^ some functions, 0o is the uniform medium order param- 
eter. A dimensionless variable x = x/L has been introduced with L being the characteristic 
length scale so that dx^{x) ~ 1. The terms with two and four spatial derivatives are then 

'm|<l>|y/^ [ l-2(7-ln2) 



27r 





|V^<I>| 


2 


8 m ( 




l<l>l 


|2£ 



1 + '-e], (16) 



(|V,$ 


2^2 




l<l>l' 


4 



+ - (17) 



where Z41, Z22 are unknown functions. We may determine magnitude of the Z's by con- 
sidering a weakly inhomogeneous configuration, $(x) = 0o + S(f){x) with \6(f)\ ^ 0o- On one 
hand, in this case (fT7|) reduces to 



where terms of higher order in 50 are omitted. On the other hand, 50 derivative terms may 
be obtained by expanding (jl]) with J = and (j){k) = 0o + 5(j){k) through the 4th order 
in 6(j){k) and in the transfer momentum, k. This produces the following operators in the 
effective Lagrangian 

m0o \ (^^ |W0(a;)p ^ (|W0(x)n^ ^ \^'5<P{x)f , .^g^ 



6^{k) 




FIG. 1: Diagrams that produce the LO contributions to the operators in (I19p . The fermion 
propagators in the </>o background are the sohd hues, the 5(f){k) insertions are the dashed hnes. The 
arrows indicate the direction of particle number flow. Only normal contributions are shown. 

where C22, C41 and C12 are some coefficients. Shown in Fig. ([1]) are the diagrams that give 
LO contributions to f|T9|) . The coefficient at k"^ \5(f){k)\'^ is given by a divergent integral 
which in the dimensional regularization is 0{e^^), while /c^|(50(/c)p and /c'^ |50(A;)|^ terms 
are multiplied by convergent integrals and, thus, have coefficients of order one. So, we see 
that C22, C41 are of order one, while C12 is 0{e~^). 

Ignoring an unlikely possibility that a Z{1) and a Z ($/0o) 1 ^ ^ 'Po niay be of different 
order in e and combining this with an observation that the sum of Z^i and Z22 is given by 
the 0{1) coefficient of the {A^Y term of Vefr($, ^o? ^)) which was already mentioned in the 
end of the previous section, we conclude that Z41 and Z22 are of order one. 

More generally, an operator A;" |(50(A;)|' with n > 2 and / > 2 will have a coefficient 
given by a convergent integral 0{1). Then all the higher order derivative terms will have 
coefficients 0(1). 

We observe presence of two length scales in the derivative terms: 

~ {em<P,)~'/^ ~ (m/x)-V2 ^ e-''\ (20) 
which is proportional to the inter particle separation, and 

~ (m0o)"'/'~l, (21) 

which is superfluidity scale (coherence length). 

Note that L^/Lgf ~ e^^^"^ ^ 1 in the e expansion. In the case L > the four and 
higher order derivative terms are suppressed by at least two powers of e. 
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So, for the situations where one is interested in large enough scale phenomena so that 
L > {m jj)'^^"^ and typical time scale T > yU^^ the effective Lagrangian of unitary Fermi gas 
to NLO in the e expansion is 



Le//[$(a;),/i] 



m |$| 
27r 



d/2 



Am 
|$| 



|V<I> 



1 /i 

2|$|2e ^ 7 
7-3(7 + ln2 



1 - 2(7 - In 2) 
1 + — -e 



4 



m |$| 
27r 



d/2 



6 



(22) 



where 7 ~ 0.57722, is the Euler-Mascheroni constant and C ~ 0.14424. 

Now let us quote the NLO effective Lagrangian for the polarized unitary Fermi gas 7^ 
valid for L > (m/i)"^/^ and T > (the effective potential was derived in 13|, [l5|). The 
derivative terms are deduced by the same line of argument as in the symmetric case. 



Le//[$,/i,5/i] 



X 



1 + 



4m 

1 -2(7 -In 2) 



IV$| 



i /i 

2\We ^ 1 



m |$| 
27r 



d/2 



e + |e(/i - 1) (h - \n{h + VW^] 

i + I^^(l±M,^,i_e(ft-i))-3C. 



m |$| 

m |$| 

m |$| 
2n 

m |$| 
2n 



d/2 



6 



d/2 



e{h-i) 



d/2 



e(/i 



6\<t>\ 

~2A 
6|$| 



{h - 1)2(/^ + 2) ln(/.2 - 1) + ^ (1 + A(i3 _ 6 7) 



(l + ^(ll-67))-^ln(/i + l 



d/2 



e(/.-i).|$|(^(<i.) + 7^($)), h=^^ 



(23) 



dx 



dx 



dy 



[/(x)-x][/(t/)-t/] 

/(a:) /(y) 



bXa^^y) - Vj(a;,l/)^ -xy]. 



/(x) + a; /(y) - ?/ /— 

f[x)f{y) 



where 
^($) 

7^($) 

j(a;,2/) 
A 

Let us check if the LO terms of the superfiuid mode Lagrangian may be obtained from 
(122]) . Setting $(x, t) = p{x, t) exp 2 i /3(x, t) and using the solution to the equation of motion 
for the radial mode p = p(/5; ^) we get 



fix) + f{y) + {x + y)/2, k{x, y) = f{x) - f{y) - {x + y)/2, f{x) = V^^, 
{6pmf - 1. (24) 



n|V/?(f)p 
2m 



(25) 



where 



n 



2tt 



-) 



1 + 



4 



7 



+ 2 Log 2 



with C = 0.14424, is the NLO number densit 
terms of the superfluid mode Lagrangian 17 



/io = Ai/e, 



(26) 



This reproduces previously found LO 



Let us comment on the relation of to the interparticle separation. Using (!26l) we 



observe that ~ 



int 



n 



'^1'^ is the interparticle separation. Numer- 



ically ~ Lint except for very small e which is not the region we are ultimately interested 
in, anyway. For example, for e = 0.2 ^i/^+i/f^ ~ 0.3 to NLO, while the four derivative term 
is suppressed by = 0.04. 



III. UNITARY FERMI GAS VORTEX STRUCTURE 

The U{1) particle number symmetry is spontaneously broken in the superfluid ground 
state of unitary Fermi gas and. therefore, stable vortex corrflgurat.ons are expected and have, 
in fact, been observed [5]. So, as a first application of the functional obtained (and as a 
sensibility check of the formalism) we will investigate the structure of a superfluid vortex in 
the symmetric (5/i = 0) unitary Fermi gas. The problem has been addressed in, e.g., [3] 
using a different method. 

Using Lagrangian fl22l) consider a single vortex configuration of unit winding number 
and set $(x) = p(r)e*^ with x = {r,6,...}. We will assume that a vortex filament is two 
dimensional in any number of space dimensions. The resulting equation of motion has been 
solved order by order in e to NLO. The LO, NLO contributions and their sum with e = 1, 
the prediction for 3 dimensions, are shown in Fig. (1). In the region near the vortex core 
for r ~ (uKpo)'^^'^ <^ Lfj_, the derivative expansion will break down. However, since in a 
vortex $(r) — » as r — 0, the error in the contribution to, for example, free energy from 
the region r < (m0o)~^^^ is insignificant. The typical size of a single vortex configuration. 
To, defined by p(ro) = p(r = oo)/2 is tq = .43 (myu)^-*^/^ at LO and ro = .45 (mp)"-'^/^ at 
NLO. In order to relate these results to experimentally relevant quantities, ( l22l) was used to 
study $(2;) profile in a harmonic trap where 

^ (1 - ^) • (27) 
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5 10 15 

FIG. 2: The single vortex profile obtained from Lagrangian ()22|) . The LO piece is in green, the 
NLO is in blue, their sum with e = 1 is the red curve. The NLO (e = 1) background value of the 
order parameter is shown for reference (/io = Ai/e)- 

The trap sets the scale, so, even the LO derivative term is suppressed by L'j^/ L"^ ~ A^~^/^, 
where N is number of particles in the trap. Derivative term becomes relevant only near 
the boundary of the cloud where L ~ (m (/i - U(f)))"^/2 [17]. This illustrates validity 
of Thomas- Fermi approximation, a well known fact. Then, using Thomas- Fermi solution, 
$(r) = 00 (yU [1 — r^/L^]) , where 0o(/^) is given in ([9]), one expresses the size of a vortex in 
terms of the trap size and particle number. In the spherical trap case we get ~ 0.20(1 + 
eO.27) L {N e)~^^^. In the cylindrically symmetric case tq ~ 0.23(1+£0.21) (eV/iV)^/'^ where 
V is the total volume of the trap. The predictions for d = 3 are then 

ro = 0.25LA^~^/^ 

ro = 0.28 (V/A^)^/^ (28) 

for the spherical and cylindrical traps, respectively, where V is the total volume of the 
cylindrical trap. The two configurations are the limiting cases which span shapes of the 
trapping potentials realized in the experiments. 
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V(2 7r)^/V/io(/iom)'^/2 
-1 . 6r 

-1.7- 




FIG. 3: The NLO effective potential of the bulk medium V at (^/i ~ as a function of the order 
parameter (p {fiQ = fi/e). 

IV. NORMAL-SUPERFLUID INTERFACE IN THE POLARIZED UNITARY 
FERMI GAS 

In general, derivative terms are relevant for any inhomogeneous phenomenon with typical 
length scale of order of inter-particle distance. Such a situation arises in description of the 
structure of the superfiuid-normal phase interface encountered in the trapped imbalanced 
gas (612^0) 0,0,3. 

Shown in Fig . (3l) is the NLO effective potential of the bulk medium V at the critical 
polarization [l3|, [l5|. Using effective Lagrangian for the polarized unitary Fermi gas fl23|) one 
is to solve for the order parameter spatial profile connecting superfiuid and normal minima 
of the potential. Unfortunately, due to flatness of the LO potential [13|, ll5|] the interface does 
not exist at the LO in e and the problem is partially non-perturbative. Let us determine 
the structure of the superfiuid-normal phase domain wall using Lagrangian fl25]) with e set 
to 1. Note that £e//[<I>, yU, is not smooth as indicated by the presence of the 6 function 
terms. This appears to be an artifact of the e epsilon expansion which treats the chemical 
potential as a perturbation. Shown in Fig. (jlj) is the resulting profile ^{z). It has a small 
discontinuity in its derivatives at ^{z) = 2 /io; its width is of the order of the inter-particle 
separation, so the derivative expansion used to approximate effective action ([3]) with the 
effective Lagrangian (J23l) is valid in this case. Since we used £e//[$, A*, 5/^] at e = 1, the 
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FIG. 4: The superfluid-normal phase interface from the NLO effective Lagrangian with e = 1. 

location of the superfluid minimum and the critical value of are slightly different from 
the results found in the perturbation theory Q . The surface free energy per unit area 
is then 



cr 



dz (/:e//[<^'W,/i,5/i,]-£e//[$(^ = 00),/i,5/xJ) =0.58^/^ (^) + 



+ 0(5^), /Xo = p (29) 

where ^i^z) is the solution shown in Fig. (j4]). In the case of polarized gas in a spherical trap 

/ a ~ /io(/i mfl'' [ (30) 
Js Js 2 y/mJIL 

where n = n/ {fio rriY^'^ is the (dimensionless) density on the superfluid side of the interface; 

is the (dimensionless) integral over the interface. Note that the total surface energy is in 
the units of yu(/im)'^/^ which is the units of the bulk free energy. 

In the Thomas-Fermi approximation location of the superfluid-normal phase boundary 
in a harmonic trap is determined by 

5fi = 6fi,ix) =2fioix) il-l3e) + 0(6^), (3 ~ 0.47, /io(x) = - (l - r^L'') (31) 



13l . llSl ] which means that the interface is located at the distance 



R = Li 1- 



2/io 



(32) 



13 



from the center. To calculate l/^m/i L we use Thomas- Fermi solution for the order param- 
eter of polarized gas in a spherical trap 



$(r) 



e{R-r), (33) 



where (poifi) is the homogeneous order parameter ([9]) and the superfluid core radius, 
R, has been defined in fl5^ . The result is almost independent of polarization, 6 = 
(N| - N|)/(N| + N|); for 5 = 0.53 we get 

^ [eN)-^^'^ 0A7{1 + 0.21 e) (34) 



which for (i = 3 gives 



For N = 10^ the coefficient is 0.004 which is four times the value of 



/ a ~ fiifimfn"" f -n^'H.lQN-"\ (35) 
Js Js 

m 



V. SUMMARY 

We have calculated effective Lagrangian of the unitary Fermi gas using e expansion tech- 
nique as well as derivative expansion. The leading order derivative terms have been deter- 
mined by gauging the U{1) particle number symmetry and requiring gauge invariance of the 
theory. We have argued that in some realistic situations the derivative and e expansions are 
related: for large enough scale phenomena so that the typical length scale is not smaller than 

~ (mjj,)^^^'^ the LO derivative terms are 0{1), while the higher order derivative terms 
are at least 0{e'^) and, so, are to be neglected when working to NLO. Quasi-locality of the 
effective action functional {i.e. validity of the derivative expansion) is due to the feature of 
the e expansion that fi ^ e and is, therefore, treated as a perturbation. In this approach to 
inhomogeneity description one has to deal with differential equations (equations of motion 



for the order parameter) which is simpler than the usual Bogoliubov-de Gennes method [16 1 
requiring having to deal with non-local effective action functional. 

As an application of the effective Lagrangian obtained, we have determined the structure 
of a vortex in the unitary Fermi gas, and the superfluid-normal phase interface profile, all 
to the NLO in e. Predictions for the vortex size in spherical and cylindrical traps are shown 
in ( [28|) : surface free energy results are given by Eqs. fl29f35p . 

Some other applications of the functional include 

14 



study of multiple vortex configurations and calculation of the parameters of the vortex 
lattices observed in the unitary Fermi gas; 

study of properties of strongly interacting Fermi gas in a periodic potential (optical 
lattice), a subject of recent experiments \2l\ : 



recently new Monte Carlo simulations of trapped few fermion system at unitarity have 



been performed 
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2J]. Assuming validity of continuous medium description, 
one should consider the case Ltrap ~ (m/i)^^/^. Thomas Fermi approximation is no 
longer valid, while derivative expansion based functional ( 122]) should still be valid. For 



systems with fewer particles an alternative technique such as in, for example, 
be employed. 



25| must 
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